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Abstract 

We introduce the real space correlation function of B-mode polarization of the cosmic microwave 
background (CMB) as a probe of superhorizon tensor perturbations created by inflation. By causal- 
ity, any non-inflationary mechanism for gravitational wave production after reheating, like global 
phase transitions or cosmic strings, must have vanishing correlations for angular separations greater 
than the angle subtended by the particle horizon at recombination, i.e. 9 > 2°. Since ordinary 
.B-modes are defined non-locally in terms of the Stokes parameters Q and U and therefore don't 
have to respect causality, special care is taken to define 'causal .B-modes' for the analysis. We 
compute the real space .B-mode correlation function for inflation and discuss its detectability on 
superhorizon scales where it provides an unambiguous test of inflationary gravitational waves. The 
correct identification of inflationary tensor modes is crucial since it relates directly to the energy 
scale of inflation. Wrongly associating tensor modes from causal seeds with inflation would imply an 
incorrect inference of the energy scale of inflation. We find that the superhorizon B-mode signal is 
above cosmic variance for the angular range 2° < 9 < 4° and is therefore in principle detectable. In 
practice, the signal will be challenging to measure since it requires accurately resolving the recombi- 
nation peak of the .B-mode power spectrum. However, a future CMB satellite (CMBPol), with noise 
level Ap ~ 1/xK-arcmin and sufficient resolution to efficiently correct for lensing-induced .B-modes, 
should be able to detect the signal at more than 3a if the tensor-to-scalar ratio isn't smaller than 
> r ~ 0.01. 
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1 Introduction 



Inflation [1-3] provides a convincing microscopic explanation for the observed scalar fluctuations 
in the temperature of the cosmic microwave background (CMB). The once popular causal seeds 
models for structure formation have been ruled out since they fail to reproduce the Doppler peaks 
of the CMB temperature correlations [4-7]. In addition, the superhorizon nature of the scalar 
modes produced by inflation leaves a characteristic signature in the temperature-polarization cross- 
correlation [8] that cannot be reproduced by causally-constrained theories. The dominant mechanism 
for the generation of scalar perturbations is therefore a period of accelerated expansion prior to the 
hot Big Bang. Nevertheless, the constraints coming from the observations of scalar modes (i.e. 
temperature and E'-mode polarization anisotropies and their cross-correlation), don't rule out the 
possibility that, while subdominant for scalar perturbations, causal seeds may be a significant source 
of tensor modes. These tensor modes correspond to transverse, traceless metric perturbations or 
gravitational waves, gij = 5ij + hij, where h\ = d l hij = 0. 

i3-modes of CMB polarization are often described as a 'smoking gun' signature of inflationary 
gravitational waves [9, 10]. Imagine therefore the glorious day that we detect l?-mode polarization. 
How do we prove to the skeptic that the signal really has its primordial origin in inflation? In this 
paper we study causality constraints on the I?-mode correlation function in real space. This is in 
the spirit of Ref. [8] who studied the corresponding problem for scalar modes and identified the 
superhorizon signature in the temperature-polarization cross-correlation function. Here we extend 
their treatment to tensor modes and B-mode polarization. 

Consider concretely two qualitatively different mechanisms for the production of a stochastic 
background of gravitational waves: 

• First, symmetry breaking phase transitions after the hot Big Bang (i.e. after reheating) may 
produce a spectrum of tensor modes, e.g. [11]. Furthermore, defect models, e.g. [12, 13], also 
create vector and tensor modes that source CMB polarization. In these models, causality 
constrains any physical real space correlations to vanish on superhorizon scales. In the power 
spectrum this corresponds to an approximate white noise spectrum on superhorizon scales. 

• Second, quantum fluctuations during an epoch of inflation before the hot Big Bang produce 
a spectrum of long-wavelength gravitational waves. The spectrum is scale-invariant even on 
superhorizon scales in apparent violation of causality. 

While both mechanisms are causal, causality imposes a much stronger constraint on tensors from 
phase transitions or defects than on tensors from inflation (see Figure l). 1 The causal structure of a 
(flat) Friedmann-Robertson- Walker (FRW) universe is most conveniently understood by considering 
its line element in the form 

ds 2 = a(r) 2 [_ dr 2 + dx 2] ; (1) 

where r = / dt/a(t) is conformal time which equals the particle horizon at time t. In standard Big 
Bang cosmology (without inflation) there is a singularity at r = 0. Inflation extends conformal time 
to negative values and the time r = becomes the (non-singular) reheating surface. The singularity 

1 ln the following wc will therefore loscly refer to defects and phase transitions as 'causal theories' although 
there is nothing acausal about inflation. 
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Figure 1: Causal structure of the universe. Correlations between any local variables at any two 
spacetime points vanish if their backward light cones fail to intersect on the spacelike 
hypersurface £ corresponding to the phase transition at r = r pt [14]. On the surface of 
last-scattering at r rec this corresponds to angular separations 6 > 6 C ~ 2°. Longer range 
correlations are established during inflation at negative values of conformal time, r < 0. 



is pushed to r = — oo. Recombination takes place at T rec and conformal time today is tq. Let 
the phase transition occur at a time r p t shortly after reheating. The initial perturbation variables 
induced by the phase transition are defined on a Cauchy surface £ at r p t- For causal theories, the 
unequal time correlator of the source stress energy tensor T* satisfies the following constraint 

(T° u (0,T)T p \(r,T>)} = V r>T + T>. (2) 

At recombination the particle horizon is ~ r rec , which today corresponds to an angle 0^ or ~ 1° 
on the sky (see Appendix A). If we could observe fluctuations at recombination, their angular 
correlations should therefore satisfy 

C{6) = V 6 > 9 C = 29 hor w 2° . (3) 

We recognize that CMB polarization (being generated only by scattering of CMB photons off of free 
electrons) offers the opportunity to study correlations directly at recombination (this is in contrast 
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to temperature anisotropics which receive contributions after recombination, e.g. from the integrated 
Sachs- Wolfe effect). 

In this paper we discuss the superhorizon signature of inflationary tensor modes in the real 
space B-mode correlation function. This real space treatment has the advantage that causality- 
considerations are applied most directly. In fact, in harmonic space the difference between the 
angular power spectra Cf for inflationary B-modes and causally-generated B-modes can be more 
subtle than one would naively anticipate (see Appendix B). To see this consider the scale-dependence 
of the power spectra of tensor modes: inflation, of course, famously produces a nearly scale-invariant 
spectrum even on superhorizon scales 

P h (k,r) = A t k nt - 3 , nt«0, V k<k hor = l/T. (4) 

We give causal theories 'maximal benefit of the doubt' by allowing them to be tuned to reproduce the 
inflationary spectrum perfectly on subhorizon scales. Any difference between inflation and causal 
theories will then be encoded in a difference in the tensor mode power spectrum for scales near 
or above the horizon scale. (In some sense, this is the best causal theories can possibly mimic 
inflation; in realistic models one often finds a difference even on scales smaller than the horizon.) 
On superhorizon scales we asymptotically expect a white noise spectrum for causal theories 

pjf\k, t) = const. , V k < k hor = 1/r . (5) 

The angular power spectrum Cf is a projection of Ph(k) onto the sky. Ignoring very low multipoles 
which are dominated by polarization generated by scattering during the epoch of reionization, the 
low-£ scaling of the .B-mode angular power spectrum Cf is independent of t for a scale-invariant 
input spectrum (4); see e.g. [15]. Interestingly, the same ^-scaling is obtained for the white noise 
input (5) from causal theories (see Appendix B for an explanation of this projection effect). The 
difference in the fc-scalings of inflationary tensors and causal tensors for small k is therefore not 
reflected in a difference in the low-£ scaling of Cf as one might naively expect. Any difference 
between inflationary tensors and causally-generated tensors is encoded in the shape of the angular 
power spectrum near the peak. The precise difference will be strongly model-dependent. To avoid 
any model-dependent considerations we instead study the unique signature that inflation imprints 
on the superhorizon B-mode correlation in real space. 

To study the causality of B-modes in real space we have to address an important subtlety: 
B-modes are defined non-locally in terms of the Stokes parameters Q and U. The causality of Q 
and U therefore doesn't imply causality of .B-modes. To avoid this ambiguity we introduce in §2 
.B-modes as a causal alternative to the conventional B-modes. The B-modes will be defined locally 
in terms of derivatives of the Stokes parameters Q and U and are therefore manifestly causal by 
virtue of Q and U being causal. Like the regular B-modes, B-modes are only sourced by tensor 
(and vector) perturbations and therefore provide an equally good characterization of inflationary 
gravitational waves. However, B-modes have the crucial advantage over B-modes that causality is 
studied unambiguously. 
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2 Causality and Superhorizon Correlations 



2.1 Harmonic Expansion 

We recall briefly some basic aspects of the analysis of CMB temperature and polarization anisotropies. 
More details may be found in the following reviews [16] . 

The radiation field (temperature T and Stokes parameters Q and U) is expanded in spin-0 and 
spin- 2 spherical harmonics [9, 10] 

T {n) = ^2a T ,e m Y im (n) , (6) 

lm 

P±(h) = (Q ±iU)(h) = ^a± 2 ,&n ±2Ye m (n) , (7) 

lm 

where the unit vector h denotes the direction on the sky. The spin-2 polarization field may be 
transformed into rotationally-invariant spin-0 E- and 5-modes 



E (n) = ^2a E ,e m Ye m {h) , (8) 

im 

B {h) = ^aB,Em Y im{n) , (9) 



lm 



where 



«£/m = [ a 2,lm + a~2,lm] , (10) 

0-B,lm = -^r [d2,em ~ 0—2,lm\ ■ (H) 

The rotationally-invariant angular power spectra are 

C f Y = 21+ I ^Z^ a *X,tm a Y,lrn) , (12) 
m 

or (a* x £ m a,Y,e' m ') = Cf Y Sur5 mm ', where X,Y = T, E, B. The power spectra Cf Y are the Legendre 
transforms of the real space angular correlation functions 

C XY (9) = (X(h)Y(ti)} = V Cf Y P t (h -ti), h ■ ti = cos 9 . (13) 

By statistical isotropy the correlation functions depend only on the angle 9 between the two directions 
h and h'. 



2.2 Causal 5-modes 

For causal theories (defect models and phase transitions) the real space angular correlation functions 
of the Stokes parameters vanish outside the particle horizon at recombination, 9 > 9 C = 2#hor- 
However, since E- and S-modes depend non-locally on the Stokes parameters this causality constraint 
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does not have to hold for E- and £>-mode 
the following real space quantities 



correlations. To avoid this ambiguity, it is useful to define 



E(h) = 
B(h) = 



2i 



3 3 p+ - 6 6 P- 



(14) 
(15) 



Here, 6 and 6 are the standard spin-raising and spin-lowering operators, respectively [10]. Acting 
on the spin-2 spherical harmonics they give 



6 6 +2 Y em (h) = n e Y lm {h) , 6 6 _ 2 l^ m (n) = n e Y tm {n) , 



where 



ne = 



{£ + 2)1 
C-2)!J 



1/2 



E and B therefore have the following harmonic expansions 

E(n) = J2 a E,£m Y e™(n) , 

em 

B(h) = Y, a i 

where 



Yem{n) , 



a (E,B),em = n l a {E,B),im- 



(16) 
(17) 

(18) 
(19) 

(20) 



Crucially, E and B are causal by virtue of being defined in terms of derivatives of Stokes param- 
eters (which are manifestly causal). In addition, S-modes are non-zero only for non-zero S-modes, 
so B is just as good a measure of primordial tensor modes as B. However, B has the advantage 
over B that causality is checked more directly; {B has the disadvantage that its noise spectrum is 
very blue; see §3). Finally, we note for completeness, that in the small-scale limit, 6 — > — (d x + id y ), 
6 — ► — (d x — id y ), Equations (14) and (15) become 



E(h) -► -{d 2 x - d 2 y )Q - 2d x d y U = V 2 E, 
B{h) -{dl-8 2 y )U + 2d x d y Q = V 2 B, 



(21) 
(22) 



where V 2 is the two-dimensional Laplacian defined in the plane of the sky. 
For the B-mode correlation function we find 



C B (6) 



47T 



(23) 



where 



Cf = n 2 Cf « tCf . (24) 

Similar expressions hold for C TE {9) and C E {6). The nf w £ 4 conversion factor to go from Cf to 
Cf has important consequences: small scales (large £) are weighted more strongly. This has the ad- 
vantage that superhorizon contributions from the reionization signal at very low-£ are automatically 
filtered out; however, it has the disadvantage that the noise spectrum is very blue (see §3). We have 
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to live with that negative aspect of B-modes, since causality can be studied unambiguously only for 
B-modes. One way to see this, is to note that the real space correlation functions of !?-modes and 
regular B-modes are related as follows 

C B {9) = (V 2 + 2)V 2 C B (9); (25) 

the differential operator (V 2 + 2)V 2 becomes n\ in harmonic £-space, (V 2 + 2)V 2 <-> n 2 . Equation 
(25) shows that causal B-modes (implied by causal Q and U), C B {9 > 9 C ) = 0, do not imply causal 
B-modes since superhorizon B-modes may receive contributions that are in the kernel of (V 2 + 2)V 2 . 
We therefore restrict our discussion of real space correlations and causality to l?-modes. 

2.3 Superhorizon Tensor Modes 

For a causal theory we expect 

C* ? (9) = , for 9>9 C = 29 hor , (26) 

where X, Y = T,E,B and #hor ~ 1° is the angular separation corresponding to the horizon at 
recombination. Note that (26) is a model-independent constraint that any causal theory has to 
satisfy. Temperature fluctuations T and polarization i?-modes are sourced dominantly by scalar 
(density) perturbations. The litmus test for scalar superhorizon correlations therefore is the TE 
cross-correlation on scales corresponding to 9 > 2° 

C TB (9)^0, for 9>9 C ^2°. (27) 
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This is related to but not the same as the signal discussed in [8]. 2 Here we propose the analogous 
test for tensor superhorizon correlations 

C S (0)^0, for 9>9 C ^2°. (28) 

Figure 3 shows the superhorizon signal from inflation. Its characteristics are a positive peak at 
9 > 2.5° and a negative peak at 2.0 < 9 < 2.5°. 
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Figure 3: C B {9) = (BB}(9): real space correlation function on superhorizon scales. C B {9 > 2°) ^ 
is a unique signature of inflationary tensor modes. 



2 Ref. [8] phrased the causality constraints on scalar modes directly in terms of correlations of the Stokes pa- 
rameters Q and U and their cross-correlations with the temperature fluctuations: for causal theories C®Q(8), 
C uu (9) and C T ®(9) vanish for 9 > 2°. The proof of the superhorizon nature of adiabatic scalar fluctuations 
is now often associated with the negative peak of the TE angular power spectrum Cj E for 50 < i < 250 
(e.g. [17, 18]), since causal theories (or inflationary models with a significant isocurvature component) tend 
to predict positive correlations for those multipoles. However, in light of the present discussion superhorizon 
TE correlations cannot always be unambiguously identified. For TE correlations this problem is absent. 
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3 Detectability of the Signal 



In §2.2 we introduced .B-modes as an alternative to the conventional .B-modes as a tensor mode 
diagnostic. A treatment in terms of .B-modes is required in order to unambiguously discuss super- 
horizon correlations in real space. However, when quantifying the detectability of the inflationary 
superhorizon signal we have to face the major disadvantage of the .B-modes: their noise spectrum 
is very blue. 

3.1 Noise Spectra and Smoothing 

The noise in the Stokes parameters and in the .B-modes derived from them is usually assumed to 
have a white noise power spectrum (multiplied by the Gaussian window function corresponding to 
the beam of the experiment) 

N t B = A%^^H, (29) 

where Ap ]C g is the constant noise per multipole and = 0.425 #fwhm is the size of the beam. 
Reference sensitivities for representative future CMB polarization experiments are given in Table 1. 
We include in the effective noise spectrum a B-mode contribution from the gravitational lensing 
of B-modes. Lensing B-modes have an approximate white noise spectrum (for £ < 700) and can 
therefore be combined with the instrumental noise [19]. Without lensing-cleaning (delensing) the 
extra noise associated with lensing .B-modes is Ap ; i ens i ng ~ 5 /iK-arcmin. As we will see below, 
this 'lensing noise' would dominate over the inflationary signal, so delensing of the measured B- 
modes is an important constraint on any experiment proposed to measure the superhorizon .B-mode 
signal. Representative values for the residual noise Ap e ff after delensing are shown in Table 1. 
The precise values will depend somewhat on the resolution, en-,, and the instrumental noise, Ap, 
of the experiment. Further details on residual noise after delensing may be found in Ref. [19, 20]. 
The effective noise should also include residual noise from astrophysical foregrounds like polarized 
synchrotron and dust radiation from our galaxy [21]. We imagine that this is included in Apefj; a 
more accurate treatment of residual noise from foregrounds and the effects of foregrounds on the 
delensing is beyond the scope of this paper. 

In harmonic space, B- modes and B- modes are related by the conversion factor ri£, see Equation 
(17). White noise in the Stokes parameters (or B- modes) therefore becomes colored noise in the 
.B-modes 

N? = 4 A P , cff e^l ~ t A% cS e^K . (30) 

Because of the n|— factor in the relation between .B-modes and -B-modes, the noise power spectrum 
for .B-modes diverges as £ 4 . The steep blueness of the noise spectrum means that the noise level is 
dominated by small-scale noise. To regulate the divergence in the noise we therefore smooth both 
the signal and the noise with a Gaussian smoothing function of width a s . The smoothed spectra are 

Cf = C-B e -/(/+l)of j ftB = N B e -W+l)al (31) 

Since a s » £>"b, we may in practice ignore the finite resolution of the beam. We also define the total 
signal as 

Cf^Cf + N*. (32) 
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Reference Exp. 


Ap [//K-arcmin 




Ap jC ff [/iK-arcmin 




J sky 


Planck 


10.00 




0.7 


Balloon 


6.00 




7.5 




0.1 




3.00 




4.5 




0.1 


CMBPol 


2.00 




3.0 




0.7 




1.00 




1.7 




0.7 




0.50 




1.0 




0.7 




0.25 




0.7 




0.7 


no noise 








0.7 



Table 1: Reference sensitivities for future CMB polarization experiments. The instrumental noise 
level is denoted by Ap. Shown are also representative values for the effective noise levels 
after delensing Ap e g- [19, 20]. 

In Figure 4 we show the dependence of the real space superhorizon correlations on the smoothing scale 
i s = l/cr s . Although smoothing decreases the signal it increases the signal-to-noise (see Table 2). 
In Figure 5 we show the sensitivity to the noise level for a fixed smoothing scale. We also show the 
cosmic- variance limit. Since the signal is much larger than the cosmic- variance, there is in principle 
no obstacle to measuring it given sufficient instrumental sensitivity. However, in practice, the strong 
dependence of the signal-to- noise ratio on Ap e g and £ s is the major limitation to measuring the 
superhorizon i?-signal from inflation. 



C^{6) [mK 2 ] 




e 



Figure 4: Smoothed signals in real space. The smoothing scales shown range from £ s = 400 (blue) 
to 4 = 800 (red). 
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Figure 5: Signal and noise in real space (for a fixed smoothing scale £ s = 400). The noise levels 
shown are Ap e g = 0, 1, 2, 3 /iK-arcmin. Since the signal is much larger than the cosmic- 
variance, there is in principle no obstacle to measuring it given sufficient instrumental 
sensitivity. In practice, residual noise after delensing provides a significant constraint. 



3.2 Superhorizon Signal 



We now quantify the expected signal-to-noise for future experiments. We will focus specifically on 
the capabilities of a future CMB polarization satellite (CMBPol) [22, 23]. 

We bin the superhorizon signal into angular intervals b = {9i =t ^A#}, where 9{ > 9 C and define 
a signal vector S with components (S), = C B {9i). In practice, we use iVbin = 10 bins between 
#min = 2° and # max = 4°. The different angular bins are of course correlated. The covariance matrix 
C has components 

2 



2£ + 1 
4tt 



cov 



(Cf) 2 Pe(cos9 i )P e {cos9 J ) 



(33) 



where in harmonic space the covariance matrix is diagonal if the effects of partial sky coverage are 
ignored 



cov 



pB pB 



B\2 



(2e + i)W~ e > 

We define the following measure of the signal-to-noise on superhorizon scales 

sc- 1 s= ^(smc- 1 )^),. 



s_ 

N 



(34) 



(35) 



i,j e b 



In Table 2 we evaluate the signal-to-noise for different smoothing scales £ s and effective noise 

1/2 / r \ 1/2 

levels (instrument noise plus lensing residuals) Ap j(3 g = 1 /xK-arcmin x (^y) X I ^j- I . For 
illustration, we assume a fiducial tensor-to-scalar ratio of r = 0.1 (as expected for m 2 4> 2 inflation 
and consistent with current upper bounds [24]) and a sky fraction f s ^ y = 0.7. Table 2 may then be 
read as the result for Ap e g = 1 //K-arcmin. However, we note that 



S_ 

N 



/sky F\X, £ s ) , 



where x 



A 



P,cff 



(36) 
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500 


600 


700 


800 


900 


1000 


5 
N 


31.6 


27.3 


13.2 


4.2 


1.8 


< 1 


< 1 



'jM 1 / 2 v ( W 
vO.l/ A \ 0.7 



1/2 



Table 2: Signal-to- noise for a CMBPoI-like experiment, Ap e g = 1 /uK-arcmin x 

as a function of smoothing scale £ s for a tensor-to-scalar ratio r, sky fraction / s k y and 
number of bins iVbin = 10. 



i.e. the signal-to-noise ratio only depends on the combination x = ^ s . Hence, the results in 
Table 2 can be translated to different values of r by appropriately rescaling the noise level Ap e g. 
The function F(x) is shown in Figure 6 for different values of £ s . The signal-to-noise ratio for a 
no-noise cosmic-variance limited experiment, S/N ~ 50, is insensitive to the smoothing. 



S_ 
N 




Figure 6: Signal-to-noise as a function of total effective noise Ap e g (instrument noise plus lensing 
residuals) and tensor amplitude r. Shown are the results for three different smoothing 
scales. The dashed line indicates S/N = 3. 

Finally, in Figure 7 we illustrate the prospects for 'signal detection', S/N > 10 (~ 10<r) or 
S/N > 3 (~ 3a), as a function of Ap e g and r: 

(a) if r = 0.1 (0.01), then the experiment has to allow Ap je g < 3.2 (1.0) /iK-arcmin for S/N > 10, 

(b) if r = 0.1 (0.01), then the experiment has to allow Ap i6 g < 6.5 (2.1) /uK-arcmin for S/N > 3. 

We see that even for a relatively large tensor-to-scalar ratio, a detection of the signal will require 
a satellite mission like CMBPol. For smaller r, the signal-to-noise decreases roughly as S/N oc r (for 
fixed noise Ap e ff) and becomes proportionately harder to detect. However, if r > 0.01 (the target 
range for a future CMBPol mission [22]), then CMBPol would detect the real space super horizon 
signal at more than 3a if Ap e ff ~ 2/xK-arcmin (corresponding to instrumental noise Ap ~ lfiK- 
arcmin; see Table 1). 
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Figure 7: 'Signal detection' (S/N > 10 or S/N > 3) as a function of total effective noise Ap je g 
(instrument noise plus lensing residuals) and tensor amplitude r: (a) if r = 0.1 (0.01), 
then the experiment has to allow Ap e g < 3.2 (1.0) /uK-arcmin for S/N > 10 (~ lOc); 
(b) if r = 0.1 (0.01), then the experiment has to allow Ap e g < 6.5 (2.1) /iK-arcmin for 
S/N > 3 (~ 3a). 

4 Conclusions 

In this brief note we pointed out that 5-mode correlations on angular scales > 2° are an un- 
ambiguous signature of inflationary tensor modes. Since ordinary B-modes are defined non-locally 
in terms of the Stokes parameters and therefore don't have to respect causality, special care had 
to be taken to define causal I?-modes for this analysis. The signal can in principle be used to dif- 
ferentiate between S-modes generated by inflation and I?-modes arising from causally-constrained 
phase transitions or cosmic strings. The unambiguous identification of inflationary tensor modes is 
crucial since it relates directly to the energy scale of inflation. Wrongly associating tensor modes 
from causal seeds with inflation would imply an incorrect inference of the energy scale of inflation. 
In practice, we found that it will be challenging to measure the superhorizon l?-mode signal since 
it requires accurately resolving the recombination peak of the i?-mode power spectrum. However, 
with a future CMB satellite (CMBPol) the signal should be detectable if the tensor-to-scalar ratio 
isn't too small. 

We should emphasize that our conclusions are extremely conservative since our analysis assumed 
that only the superhorizon scales can be used to distinguish the inflationary signal from the signal 
associated with causal seeds. In explicit physical models for 5-modes from cosmic strings or phase 
transitions one never finds this situation, but there is always a significant difference also in the 
subhorizon signal, e.g. [11-13]. For example, while the inflationary B-mode signal peaks on the 
horizon scale, l?-modes from causal seeds generically peak on a slightly smaller scale. For explicit 
models this difference can be used to distinguish inflationary U-modes from 5-modes associated 
with causal theories even with Planck [25] . However, the subhorizon difference between inflationary 
£>-modes and causal i?-modes is model-dependent and therefore not universal. In this paper we 
have avoided this model-dependence by focusing exclusively on the unique superhorizon signature 
of inflationary i?-modes. 
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A Particle Horizon and Causality 

In this Appendix we compute the angle subtended by the comoving horizon at recombination. This 
is defined as the ratio of the comoving particle horizon at recombination and the comoving angular 
diameter distance from us (an observer at redshift z = 0) to recombination (z ~ 1090) 

Shor = ^ • (37) 

(La 

A fundamental quantity is the comoving distance between redshifts z\ and Z2 

^-^f^r 1 ^'- (38) 

The comoving particle horizon at recombination is 

<4or = Tree ~ Tj « l(z rcc , Oo) . (39) 

In a flat universe, the comoving angular diameter distance from us to recombination is 

d A = T - T rec = 1(0, Z rcc ) . (40) 

The angular scale of the horizon at recombination therefore is 

n _ ^hor l(z Tec , oo) r a-i\ 

#hor = —7— = r- • (41) 

d A 1(0, Z rcc ) 

Using 



H(z) = H ^Jn m (l + zf + 7 (1 + zY + n A , (42) 

where Vt m = 0.27, £l\ = 1 — Q m , J7 7 = Q m /(1 + z cq ) and z cq = 3400, we can numerically evaluate 
the integrals 1(0, z TCC ) and l(z rec , oo), to find 

0hor = 1.16° . (43) 

Causal theories have vanishing correlation functions for 

9>6 C = 20 hor = 2.3° . (44) 
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B CMB Angular Power Spectra 



To gain some intuition for the possible degeneracy between the CMB angular power spectra for 
inflationary l?-modes and for causally-generated U-modes after reheating we consider simple models 
for the power spectra of tensor modes: 

• Inflation predicts a nearly scale-invariant spectrum of tensor modes on superhorizon scales 

P h (k,n) = A t k nt - 3 , nt«0, V kn<l. (45) 

The redshifting of gravitational waves when they enter the horizon leads to the following 
scaling relations for the gravitational wave spectrum at recombination (r rec = l/k ICC ), e.g. [15] 



A h (k,T ICC ) = k 3 P h (k,T rec ) oc k nt < 



( 1 

k~ 4 
k- 2 



cq 



(46) 



^eq ^ k 



Here, r eq = l/k eq is the time of matter-radiation equality. Numerically, 



50 ko , k e q = 130 ko = 2.7 k TCC . 



(47) 



Of course, in reality is a smooth function of k. A useful fitting function for a smooth 
spectrum that asymptotically goes as k ni for k < k+ and k~ n2 for k > k* is 



A h (k, r*) = k 3 P h (k, r*) = A y— 



cosh 



-A 



y = 



k_ 



(48) 



• On the other hand, any causal mechanism to produce tensor modes (like global phase transi- 
tions or cosmic strings in the early universe) may be tuned to have a scale-invariant spectrum 
on subhorizon scales at recombination, but by causality we have no reason to believe that it 
will be scale-invariant on superhorizon scales. Instead it seems reasonable to expect a Poisson 
spectrum (corresponding to uncorrelated white noise) on superhorizon scales 



P^\k) = const. , k < e k TCC , e > 1 . 



(49) 



The fudge factor e parameterizes the peak of the causal tensor mode spectrum (in physical 
models e > 1, e.g. e = 3.7 in [11]). The spectrum at recombination therefore is 



r k 3 



A^ ^ (&, T rec ) — k 3 P^ \k,T TCC ) oc < k ^ e k rec < k < ek, 



•„3 p(c) 



k <C 6 kr-. 



cq 



(50) 



k- 



£ &eq < k 



Naively, one might think that the difference in the superhorizon power spectra Ph{k < k TCC ) would 
show up as a difference in the low-^ CMB power spectra Cf . We now show that this is not the case. 
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B.l Line-of- Sight Formalism 

In the line-of-sight formalism the angular power spectrum may be written as [26] 



Cf Y = (4tt) 2 J k 2 dk t\(k± A xe (k)A Y e(k\ (51) 



Spectrum Anisotropies 

where 



/"TO 

A xe (k) = / dr S x (k,r) P X e(k[T - r]) . (52) 

J0 N v ' S v ' 

Sources Projection 

For B-mode polarization the source function is 

S B (k,r) = -g*, (53) 

where g(r) is the visibility function. For a Gaussian visibility function and assuming slowly varying 
sources over the surface of last-scattering (of width Ar rec ) one furthermore finds 

* oc A(r rcc )AT r ec • (54) 

The signal is therefore determined by the time-derivative of the gravitational wave amplitude at 
recombination, /i(r rcc ) oc fc/i(r rcc ). The projection factor for S-modes is 

P Bt ( x ) = 2j' e (x) + ^-, (55) 

where x = k{jQ — r). After approximating the integral over conformal time r by the value of its 
integrand at recombination one finds [15] 

A Be (k) oc P B£ [k(ro ~ r rC c)]A(r rcc ) e -( feAr -) 2 / 2 , (56) 

and 

Cf oc / k 4 dk P h (k, r rcc ) P B e[k(T - r rcc )] 2 e ^ kAr ^ 2 . (57) 

B.2 Asymptotic Scale-Dependence 

The projection factor P B e(x) is peaked at x ps £ and for large k scales as P^ e ~ k~ 2 . This implies 
that the behavior at last-scattering of the mode with wavenumber k « £/{tq — r rcc ) dominates the 
contribution to Cf , unless the integrand grows faster with k than the tail of the projection function 
(this pathological case will in fact be relevant for 5-modes from causal theories). 

• For inflation, one can make the following instructive approximation 

Cf a k 5 P h (k,r rec )\ k ^ /{To _ Tiec) J dlnx P m [x] 2 . (58) 
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From (46) we may then infer that the polarization power spectrum satisfies the following 
scalings 



£ £ <C £rcc 



£{£ + l)Cf oc < 



eq 



(59) 



r 4 £ A < £ 

• We have given causal tensors 'maximum benefit of the doubt' and assumed that pff\k) is 
the same as for inflation for subhorizon scales at recombination, k > k rcc . The expected Cf 
spectrum is then the same for £ > £ rec . Let us compare the spectrum for £ <C £ re c- The 
approximation (58) is now invalid since the integrand grows faster than the decay of the tail 
of the projection function. Instead, the integral is 



Cf 



for £ < £ rec , where x rec = 
of £, we find 



/ x^dx Pbi[x] 2 + small correction , (60) 
Jo 

« 50. Since for low £ this converges to a constant independent 



+ \)Cf oc£ z £<£ T 



(61) 



just as for inflation! Hence, we conclude that although the power spectra Ph(k), (46) and 
(50), are very different for superhorizon scales, k < k TCC , the projection effects make the Cf 
spectra look very similar even at low multipoles, £ <C £ re c- 

Since we assumed that the spectrum on subhorizon scales was tuned to be very close to the 
spectrum from inflation and we just showed that on superhorizon scales projection effects make the 
scaling of the Cf spectra virtually identical, we suspect the causality constraint to be encoded in a 
difference in the spectra around £ rN -' -tree* 

The precise difference in the Cf spectra at £ 

^ £ycc is model-dependent, so instead we focused in 
this paper on a model-independent signature of superhorizon i?-modes in the real space correlation 
function. 
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